Quark matter in neutron stars within the Nambu - Jona-Lasinio model and confinement 
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The quark matter equation of state (EOS) derived from the standard Nambu - Jona-Lasinio (NJL) model is 
soft enough to render neutron stars (NS) unstable at the onset of the deconfined phase, and no pure quark matter 
can be actually present in its interior. Since this is a peculiarity of the NJL model, we have studied a modified 
NJL model with a momentum cut-off which depends on the density. This procedure, which improves the 
agreement between QCD and NJL model at large density, modifies the standard NJL equation of state, and then 
it is potentially relevant for the stability analysis of neutron stars. We show that also within this approach, the 
NS instability still persists, and that the vacuum pressure, as a signal of quark confinement, has a fundamental 
role for the NS stability. In this respect, our conclusions point to a relationship between confinement and NS 
stability. 

PACS numbers: 97.60.Jd, 21.65.+f, 12.38Aw, 12.38Mh 



I. INTRODUCTION 

In the core of astrophysical compact objects, like neu- 
tron stars (NS) or proto-neutron stars, nuclear matter is ex- 
pected to reach a density which is several times nuclear sat- 
uration density. Calculations based on microscopic Equation 
of States (EoS), which include only nucleonic degrees of free- 
dom, show that the central density of the most massive neu- 
tron stars can be from seven to ten times the nuclear satu- 
ration density In such configurations the nucleons are 
closely packed, and to consider them as separate entities be- 
comes highly questionable. Unfortunately, it is difficult to cal- 
culate accurately the transition point from nucleonic to quark 
matter, and only rough estimates have been given in the liter- 
ature 01 . The microscopic theory of the nucleonic Equation 
of States has reached a high degree of sophistication, and dif- 
ferent many-body methods have been developed. They show 
a substantial agreement among each others yfl, and therefore 
the main uncertainty of the transition point stays on the quark 
matter EoS. Assuming a first order phase transition, as sug- 
gested by lattice calculations, one can use different models 
for quark matter to estimate the transition point, and calculate 
the compact object configuration. This approach has been fol- 
lowed by several authors, and a vast literature exists on this 
subject. The applications of quark matter models to the study 
of the NS structure have used different versions of the MIT 
bag model [4j], the color-dielectric model [5] and different for- 
mulations of the NJL model |[g]. In general, it seems that the 
maximum mass of NS that contain quark matter in their inte- 
rior is bounded to be less than 1.6-1.7 solar masses. However, 
more recently it has been shown in [7] that if one corrects the 
MIT bag model by introducing additional terms suggested by 
perturbative QCD, the maximum mass can reach values close 
to 1.9 solar masses, similar to the ones obtained with nucle- 
onic degrees of freedom. Despite the similarity of the results 
on the value of the maximum NS mass, the predictions on 
the NS configurations can differ substantially from model to 
model. The most striking difference is in the NS quark matter 
content, which can be extremely large in the case of EoS re- 
lated to the MIT bag model or the color-dielectric model, but 



it is vanishingly small in the case of the original version of the 
NJL model [6, 8]. In the latter case it turns out that, as soon as 
quark matter appears at increasing NS mass, the star becomes 
unstable towards collapse to a black hole, with only the pos- 
sibility of a small central region with a mixed phase of nucle- 
onic and quark matter. This result can be quite relevant for the 
physics of NS, since the NJL model is the only model which 
is based on phenomenological low energy data, i.e. on hadron 
properties. The main drawback of the model is the absence 
of confinement, since the gap equations which determine the 
quark masses as a function of density cannot incorporate a 
confining potential. One then assumes that the chiral phase 
transition marks also the confinement transition, as indicated 
by all lattice calculations. More recently a confining potential 
has been introduced in the NJL model [9], which is simply 
switched off at the chiral phase transition. In this case indeed 
no sharp instability of the NS is observed at the onset of quark 
matter in the central core. This is suggestive of a connection 
between the presence of confinement and the possibility of NS 
with a quark matter central core. However it is not clear if the 
instability of the NS is related or not to confinement for (at 
least) two reasons. Indeed, the introduction of the confining 
potential requires several parameters and, moreover, the NJL 
model at large density suffers of another drawback pointed out 
in [ 10]. In fact, it has been shown that the standard NJL model 
is not able to reproduce the correct QCD behavior of the gap 
for large density, and therefore a different cut-off procedure at 
large momenta has been proposed. More precisely, a density 
dependent cut-off has been introduced, and this strongly mod- 
ifies the standard NJL model thermodynamics. Therefore, the 
stability analysis of NS by the new EoS, which follows from 
this modified treatment of NJL model at large density, is a 
preliminary step to understand whether confinement is an es- 
sential ingredient in the stabilization mechanism. 

It is the purpose of this paper to clarify this point and to 
identify the origin of the NS instability at the quark onset 
within the original NJL model, which sharply distinguishes 
this model from all the others. We show that with the density 
dependent cut-off procedure the NS instability still persists 
and that the vacuum pressure, as signal of quark confinement, 
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has a fundamental role for the NS stability, as yet observed in 
the MIT bag model. In this respect, our conclusions point to 
an indirect relationship between confinement and NS stability. 



II. NJL AT LARGE DENSITY 

The NJL model provides a good phenomenological descrip- 
tion of low energy QCD based on the Lagrangian (for the two 
flavor case) 



&NJL = iyffd^y + g 
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As pointed out in 11011 . the model is not able to reproduce 
the correct behavior of the solution of QCD gap equation at 
large density. Indeed, at zero temperature, the relevant dynam- 
ical degrees of freedom have momenta in the shell \p\ ~ fl + 8, 
where jj. is the chemical potential and 8 is the cut-off evalu- 
ated from the Fermi momentum. Therefore the standard cut- 
off A on the momentum, (\p\ < A), forces 8 ~ when jj. ~ A, 
and this gives an unphysical reduction of the number of dy- 
namical states of the system. 

To solve this problem a new procedure has been proposed, 
which simultaneously preserves the low energy properties of 
the theory, and allows one to write a cut-off independent gap 
equation at zero temperature and density. The idea is the intro- 
duction of a cut-off dependence in the four-fermion coupling 
constant, g, by: 

1) keeping fixed f n to its experimental value and deriving the 
constituent mass M as a function of A, M(A), from the ex- 
pression for f n : 



fn = 3M 2 
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2) determining the coupling g(A) which solves the gap equa- 
tion, at ji = and zero temperature, 



M(A) =m + 4N/N C g(A) 



d 3 p 
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for any value of A. As shown in Fig. 1 of 111 Oh . A g(A) 
smoothly decreases with A. 

Then, according to the procedure developed in 1 1 Oil - at finite 
density, a ji dependent cut-off A(ji) is introduced, which, in 
turn, implies a jj. dependent coupling constant. In particular, 
a linear dependence in A(ju) is taken, which provides a better 
agreement with high density QCD. 

On the other hand this procedure strongly changes the ther- 
modynamics and the EoS of the system at finite density with 
respect to the standard NJL approach. In fact, the functional 
form of the thermodynamical potential Q. is not modified by a 
fl dependent cut-off, and therefore the pressure is 



p = -V- l £l(A(ii),g(Li),V,T = 0) 



(4) 



where V indicates the three-dimensional volume of the system 
and Q., in the following, will be computed in the usual mean 
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FIG. 1: The cut-off is displayed as a function of the chemical po- 
tential. Several functional dependencies have been adopted, as dis- 
cussed in the text. 



fie ld ap proximation. The particular choice of A(ji) considered 
in CM: 
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where Ao ~ 600 MeV is the NJL cut-off in the standard treat- 
ment, and c = 8/fi= 0.35 (see fiolo . introduces a discontinu- 
ity in the density p = —y(dQ./d ju), and in the corresponding 
EoS. This definition of the density p takes care of the depen- 
dence of A on ji, and ensures thermodynamical consistency. 

Since there is no compelling restriction to a specific func- 
tional form of A(ju), we replace Eq.© with some smooth 
interpolations that have no discontinuous densities and grow 
with different slopes for large values of the chemical poten- 
tial. In Fig.l we plot the cut-off of Eq. (O together with three 
different ansatz: 



Ai(M) = 



II < 400 MeV 

/9Gu-400) 2 + A2 ^> 400 MeV 



(6) 
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A 3 (M) = 
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ju < 435 MeV 
-2A li > 435 MeV 
(7) 



ju < 397 MeV 
ju > 403 MeV 



(8) 



The numerical coefficients in Eqs. ([6j, (O, dD, are adjusted 
in order to keep the curves in Fig.l either close to the one 
of Eq. (0, or close to the constant standard cut-off Aq. In 
particular, for A3, in the interval 397 MeV < pL < 403 MeV, 
a continuous quadratic interpolation between the two straight 
lines reported in Eq. ([8j, is taken. We stress that the adopted 
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FIG. 2: The energy per baryon is shown as a function of the baryon 
density (upper panel), whereas the pressure as a function of the 
baryon chemical potential is displayed in the lower panel, for dif- 
ferent choices of the cut-off. See text for details. 
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FIG. 3: The complete EOS is displayed, i.e. the pressure as function 
of the baryon density. See text for details. 



smoothing procedure is not intended to propose a method al- 
ternative to the one presented in ifioll . but it has the only pur- 
pose of making this approach suitable for our application to 
NS, and avoiding spurious discontinuities in the density. 

Considering different slopes in the [i dependence of the cut- 
off is an important issue, because this implies different growth 
behaviors of the pressure, and therefore of the density, as a 
function of the baryon chemical potential jj,b, and moreover 
different dependence of the energy per baryon on the baryon 
density pg, which is crucial to determine the critical value of 
Pb in the transition from nuclear to quark matter. 



In Fig. 2 we display in the upper panel the energy per baryon 
E/A versus Pb/Po f° r the pure quark phase (po indicates the 
nuclear saturation density, po = 0.17 fm~ 3 ), and the pressure 
as a function of jig in the lower panel for the various A,-(/i). 

For completeness, we also show the results obtained using 
a constant cut-off, as in the standard NJL model (thin solid 
line), for which the parametrization of ref.|Ht] has been used. 
As we shall see in the next section, the stability analysis of 
the NS is influenced by the various behavior of E /A plotted in 
Fig. 2. 



III. NJL AND NEUTRON STAR STABILITY 

The transition to quark matter in the core of NS has been 
studied by many authors by matching the baryonic EoS, as 
determined by microscopic many -body calculations involving 
nucleons only, and the quark EoS as determined by simple 
models which simulate QCD. The crossing of the two EoS 
in the pressure vs. chemical potential plane marks the transi- 
tion point between the two phases. At the transition the two 
phases are in mechanical and chemical equilibrium. Below 
the transition point, i.e. at lower density, the nucleonic phase 
is favored, while above the point the quark matter is favored. 
This corresponds to the Maxwell construction for a first order 
phase transition, as suggested by the indications coming from 
lattice calculations 111 ill . Solving the Tolmann-Oppenheimer- 
Volkoff (TOV) equations [ 12], one can then calculate the NS 
configuration and the corresponding gravitational mass as a 
function of the central density or of the corresponding NS ra- 
dius. In a previous work |@], where the NJL model was used, 
it was found that the quark onset at the center of the NS as 
the mass increases marks an instability of the star, i.e. the NS 
collapses to a black hole at the transition point since the quark 
EoS is unable to sustain the increasing central pressure due to 
gravity. For the nucleonic sector a microscopic EoS derived 
with the Brueckner-Bethe-Goldstone (BBG) method was used 
013I1 . The uncertainty in the baryonic EoS is relatively small 
and not relevant to the present analysis. In fact, the results are 
in line with similar calculations within the NJL model, where 
other nucleonic EoS are used [14]. For simplicity here we 
restrict the calculations to the two flavor case. 

In the calculations presented in this paper, for the hadronic 
phase we have adopted a nucleonic equation of state obtained 
within the BBG approach 01511 . using the Argonne Vis two- 
body potential 01611 . supplemented by the Urbana phenomeno- 
logical three-body force 01711 , For the quark phase, we use 
the standard parametrization of the NJL model JD], and our 
proposed prescriptions of the cut-off described in Sect.II. In 
Fig. 3 we display the complete EOS, i.e. the pressure as func- 
tion of the baryon density for the cases discussed above. The 
thick solid line displays the equation of the state for the pure 
hadron phase, treated in the BBG approach. The thin mono- 
tonically increasing curves represent the quark matter branch 
of the EOS, whereas the plateaus are a consequence of the 
first order Maxwell construction. The transition density range 
shown by the plateau depends on the adopted cut-off proce- 
dure. 
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FIG. 4: The gravitational mass (in units of the solar mass Mq = 2 x 10 g) is plotted as function of the radius (left panel) and the central 
density (right panel), for different choices of the cut-off behavior. 
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FIG. 5: The pressure is plotted as function of the energy density for 
the two-flavour case. 



In Fig. 4 we report the neutron star masses as a function of 
the radius (left panel), and of the central density (right panel) 
for the different choices of the density dependent cut-off dis- 
cussed above. For comparison also the results for the standard 
NJL model (three flavor) of ref. [6] are reported (full circles). 
In all cases one can see that at the maximum mass the plot 
is characterized by a cusp, which corresponds to the instabil- 
ity mentioned above. The plateau which appears in each of 
the curves reported in the right panel is a consequence of the 
Maxwell construction. The width of the plateau is the jump of 
the central density at the onset of quark matter. 

It is likely that if a mixed phase, according to the Glenden- 



ning construction lfl8ll . had been used, then the plot would be 
smooth, but in any case no pure quark matter is expected to be 
allowed in a stable NS. 

Therefore the improvement of the NJL model by the den- 
sity dependent cut-off does not solve the NS stability problem, 
and the origin of the instability should be related to the other 
missing dynamical ingredient, i.e. the quark confinement. 

To have an indication on the role of confinement, let us con- 
sider the pressure as a function of energy density for the NJL 
model, as reported in Fig. 5. One can see that in the relevant 
range, since the chiral symmetry is restored, the model be- 
haves like the MIT bag model (i.e. a free Fermi gas) with a 
bag constant Bnjl close to about 140 MeV fm~ 3 . It is well 
known that at such value of the bag constant the maximum 
mass within the MIT bag model is much smaller than the value 
at the maximum of the curves in Fig. 3. In general, the maxi- 
mum mass of NS increases at decreasing value of B approxi- 
mately as B- 1 / 2 fill. 

It looks likely that the same mechanism is acting in the 
three flavor case. In fact we have seen that the new cut-off 
prescription A = A(jx) tends to soften the quark EOS with 
respect to the standard NJL. Since the latter displays already 
the NS instability discussed above even with three quark fla- 
vors Jfl, we can expect an enhancement of the instability in 
the present case. Within the NJL model the value of the ef- 
fective bag constant Bnjl is dictated by the low energy phe- 
nomenology and it cannot be tuned. Indeed, the physical con- 
tent of the model demands that the pressure in vacuum is zero, 
since there is no confinement, and then the constant pressure 
Bnjl above the chiral transition is necessarily present and de- 
termined uniquely by the parameters of the model |[20ll . If 
one adds by hand a confining potential which is switched off 
at the chiral phase transition |9|], then of course the instabil- 
ity can be removed since the effective bag constant would be 
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correspondingly reduced. 

In any case, this shows that the instability is closely linked 



to the lack of confinement in the original NJL model. 



[1] M. Baldo, I. Bombaci, and G. F. Burgio, Astronomy and Astro- 
physics 328, 274 (1997). 
[2] M. Baldo, P. Castorina and D. Zappala', Nucl. Phys. A 743, 44 

(2004) . 

[3] M. Baldo, "Nuclear matter Equation of State at high density", 
Proceedings of the 5th Conference on Hadron Physics at ICTP, 
2006, Nucl. Phys. A, C. Ciofi degli Atti and D. Treleani Eds, to 
be published. 

[4] G. F. Burgio, M. Baldo, P. K. Sahu, A. B. Santra, and H.-J. 

Schulze, Phys. Lett. B 526, 19 (2002); G. F. Burgio, M. Baldo, 

P. K. Sahu, and H.-J. Schulze, Phys. Rev. C 66, 025802 (2002), 

and references therein. 
[5] C. Maieron, M. Baldo, G. F. Burgio, and H.-J. Schulze, Phys. 

Rev. D 70, 043010 (2004), and references therein. 
[6] M. Baldo, M. Buballa, G. F. Burgio, F. Neumann, M. Oertel, 

and H.-J. Schulze, Phys. Lett. B 562, 153 (2003), and references 

therein. 

[7] M. Alford, M. Brady, M. Paris and S. Reddy, ApJ 629, 969 

(2005) . 

[8] K. Schertler, S. Leupold and J. Schaffner-Bielich, Phys. Rev. C 

60, 025801 (1999). 
[9] S. Lawley, W. Bentz and A.W. Thomas, J. Phys. G 32, 667 



(2006). 

[10] R. Casalbuoni, R. Gatto, G. Nardulli and M. Ruggieri, Phys. 

Rev. D 68, 034024 (2003). 
[11] Z. Fodor and S. D. Katz, JHEP 0404, 050 (2004). 
[12] S.L. Shapiro and SA. Teukolsky, Black Holes, White Dwarfs 

and Neutron Stars (Jhon Wiley and Sons, New York, 1983). 
[13] M. Baldo, G. F. Burgio, and H.-J. Schulze, Phys. Rev. C 61, 

055801 (2000). 

[14] D. P. Menezes and C. Providencia, Phys. Rev. C 68, 035804 
(2003). 

[15] X. R. Zhou, G. F. Burgio, U. Lombardo, H.-J. Schulze and W. 

Zuo, Phys. Rev. C 69, 018801 (2004). 
[16] R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, Phys. Rev. C 

51, 38 (1995). 

[17] J. Carlson, V. R. Pandharipande, and R. B. Wiringa, Nucl. Phys. 

A401, 59 (1983); R. Schiavilla, V. R. Pandharipande, and R. B. 

Wiringa, Nucl. Phys. A449, 219 (1986). 
[18] N. K. Glendenning, Phys. Rev. D 46, 1274 (1992). 
[19] E. Witten, Phys. Rev. D 30, 272 (1984). 
[20] M. Buballa, Phys. Rep. 407, 205 (2005). 



